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X-ray shots of Cyg X-l in different energy bands and spectral states have been studied 
with PCA/RXTE observations. The detailed shot structure is obtained by superposing 
many shots with one millisecond time bin through aligning their peaks with an improved 
■ algorithm. In general, the shots are composed of a slow rise and fast decay. The shot 

structures in the different states are different. The duration of shot in the high state is 
shorter than that in the low and transition states. The shot profile in the high energy 
band is more asymmetric and narrower than that in the low energy band. The average 
hardness of shot is lower than that of steady emission in the transition and low states but 
higher than that in the high state. The time lags between the shots in higher and lower 
energy bands have been found in the different states. In transition states, the time lag is 
the largest among the different states of Cyg X-l, and it is the smallest in the low state. 
q The implications of the observed shot features for shot models are discussed. 

Subject headings: accretion: accretion disk - stars: individual (Cygnus X-l) - X-rays: 
stars 
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1 Introduction 



The black hole binary Cyg X-l has complex X-ray emission variability at all the time scales 
from milliseconds to years (cf. review by Tanaka & Lewin 1995, van der Klis 1995). On the 
short time scale of subsecond the hard X-rays from Cyg X-l show large- amplitude chaotic 
fluctuations. The studies of temporal properties of the rapid variability in the hard x-rays 
should be a hopeful method to understand the properties of the compact object. In the low 
state (LS), the variability of X-rays could be described approximately as the superposition 
of individual shots (Terrel 1972; Oda 1977; Sutherland et al. 1978; Miyamoto et al. 1988; 
Miyamoto et al. 1992). The properties of the rapid fluctuation or shot had been studied from 
ACF (Autocorrelation Function) and PDS (Power Density Spectrum) based on the supposition 
of randomly occurring shots (Nolan et al. 1981; Lochner et al. 1991). Negoro, Miyamoto and 
Kitamoto (1994) (hereafter NMK) obtained average shot profiles from a Ginga observation of 
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Cyg X-l in the low state by superposing many shots through aligning their peaks and found 
that the shot profile was approximately represented by a sum of two exponentials with time 
constants of about 0.1 and 1 s and that the average hardness during the shot is lower than that 
of steady emission. 

Quite a few models have been proposed as the mechanism for creating the rapid fluctuations, 
such as the hot spot model in the inner part of disk (Lightman & Eardley 1974; Shibazaki & 
Hoshi 1975; Rothschild et al. 1974; Miyamoto & Kitamoto 1989), the disturbance propagation 
model (Manmoto et al. 1996) and the model of magnetic flare (Galeev et al. 1979; Pudritz & 
Fahlman 1982; Nolan et al. 1981; Lochner et al. 1991) . The rapid hardening of the spectra at 
the peak obtained from the study of average shot features by NMK is difficult for the magnetic 
flare model to explain. Recently, the observed phase lags and variability coherence between 
high- and low-energy bands have ruled out the magnetic flare model and the models of hot spot, 
and it seems that the simple propagation model could qualitatively represent the properties of 
x-ray fluctuation (Nowak et al. 1997; Vauhgan & Nowak 1997). 

Up to now, the shot properties have only been studied with observations in the low state 
of Cyg X-l. The RXTE mission (Bradt et al. 1993) provided a unique opportunity for shot 
property studies in Cyg X-l with microsecond resolution observations when Cyg X-l experi- 
enced transitions from the low state to the high state, and then from the high state to the low 
state ( Belloni et al. 1996; Cui et al. 1997a, 1997b, Zhang Cui et al. 1997). In this paper, we 
present our analysis of the structures of X-ray shots in the different states of Cyg X-l with an 
improved method from RXTE observations. 

2 Observations and Data Analysis 
2.1 Observations 

There are fifteen public observations with RXTE from May 22 to August 15, 1996, which 
covered transition states (IS) and high state (HS) of Cyg X-l (Belloni et al. 1996; Cui et al. 
1997a, 1997b). The PCA observations were carried out consistently with three modes. Two 
of the modes were selected for analyzing shot structure with one millisecond time bin in this 
work: one is the Event-mode with sixteen energy bands covering 13 — 60 keV with the time 
resolution of 16 fis; the other is the Single-Bit mode including 2 — 6 keV and 6 — 13 keV energy 
bands with the time bin of 125 fis. The data from two of PCA/ RXTE observations on 1996 
Dec. 17 when Cyg X-l stayed in the low state, were also used in this study. In those data, the 
energy bands of Single-Bit modes of PCA are 2 — 5 keV and 5 — 13 keV with the time bin of 
125 /is. 

The chosen data have been further binned with a 2~ 10 s time resolution. The collimator 
response correction and barycentric correction are applied. In the study of the average features 
of shots from a bright source, such as Cyg X-l, the background can be negligible. The effect of 
dead time can be also neglected because the real count rate is not very high and the time bin 
is much larger than the 10/zs dead time (Zhang 1995). 
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2.2 Analysis Methods 



According to NMK, shots were selected by the criteria that their peak count of 1.2 — 58.4 keV 
X-rays should be 2 — 3 times larger than that of the local mean counts obtained in the 32 s 
interval and should be the maximum within the neighboring 8 s on both sides. The average 
shot profile with the 15.6 ms time bin in the 1.2 — 58.4 keV energy band and that in the five 
subdivision bands (1.2-7.3, 7.3-14.6, 14.6-21.9, 21.9-36.5 and 36.5-58.4 keV) were obtained 
separately by superposing the count profiles in each band of selected shots through aligning 
their peak bins determined by the total light curve in 1.2 — 58.4 keV. Since the differences 
between the times of shot peaks in different energy bands were not taken into account and 
that the shots near which there were other smaller shots within 8 s were not excluded in the 
NMK's procedure, the average shot structure may deviate from the truth. Li and Fenimore 
(1996) (hereafter LF) proposed a peak-finding algorithm to analyze the 7-ray burst temporal 
structure. In this algorithm, a candidate peak is selected as the bin having more count (C p ) 
than the neighboring bins (both sides), and then both sides of each candidate peak is searched 
for the bins with count C\ (at t\ < t p ) and C 2 (at t 2 > t P ) so that the condition 



is satisfied. When both C\ and C 2 are found the candidate peak is selected as the true peak 
of a pulse. This algorithm is sensitive to selecting the shots but the statistical fluctuation of 
count has more influence on its result. 

In order to analyze light curves with a time resolution down to 1 ms and low counts in each 
bin, we combined the two methods mentioned above and made necessary improvements in our 
work. To suppress the effects of statistical fluctuation, the larger time bin of 10 x 2- 10 s (~ 10 
ms) light curves is used initially to find the shot peak according to the LF method. Then each 
selected peak bin and its neighboring bins (both sides) are divided into thirty bins with time 
bin of 2~ 10 s (~ 1 ms). A shot peak is selected from the thirty bins by using NMK method 
in which the average count is used as the count of the steady emission in the interval of 1 s 
around the peak. This selecting process has been carried out separately in each energy band 
of three subdivisions between 2-60 keV. A shot peak candidate is selected as a true peak if 
it coincidences in all the three bands within 30 ms. To reduce the effects of shot overlap, we 
excluded any shots near which other shot was found within one second. Because there are time 
delays between X-ray intensity variations at different energies, the shot peaks were determined 
and aligned to obtain the average shot profiles separately in each energy band. 

The improved method used in this work, the NMK method and the LF method have been 
compared by simulations with a simulated light curve which has a similar count rate of steady 
and the shot peak as observed by PC A/ RXTE and a shot profile of exponential decay. The 
obtained average shots are displayed in Figure 1. It shows that the average shot feature obtained 
by our method is more similar to the given one than that of the NMK and LF methods. The 
average shot was not significantly changed as the value of parameter a in the selection condition 
(1) varied from 2.0 to 4.0. It means that the shots selected under the condition of a=2.0 are 
most likely to be rapid flares with time scale of 10ms, which significantly exceed the steady 
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Figure 1: Comparison of the present method and other methods, (a) Original shot, (b) 
The average superposed shot obtained by the method of Negoro et al (1994). (c) The average 
superposed shot obtained by the method of Li and Fenimore (1996). (d) The average superposed 
shot obtained by the present method. 



emission, rather than statistical fluctuations. In the following analysis, a is set to be 2.0. It 
can be seen that the obtained average shots tend to be expanded because of the selecting errors 
of peak bins due to their count fluctuation. The peak bin is excluded in our studies because 
its count has a systematic error caused by the criterion of shot peak selecting. 

We divided the selected shots into two groups with high and low peak counts, and found 
that the average shot profiles were very similar between them, as shown in Figure 2. It means 
that the shot profile in one state does not correlate significantly with the peak counts. A similar 
result was also obtained by NMK from a Ginga observation when Cyg X-l stayed in the low 
state. 

We used the full width of half maximum (FWHM) of autocorrelation to describe the duration 
of shot. The shots were divided into two groups with duration being longer or shorter than 130 
ms. The average shots of the different groups are shown in Figure 3. It can be seen that, except 
the duration, the shot features are quite similar between the long and short duration groups. 
The higher the energy bands are, the more similar the two groups are. The difference between 
the two groups can be attributed mainly to the slow variation. The rapid decay processes are 
almost the same in the two groups. The hardness evolution during the shots in the two groups 
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Figure 2: The average profiles of shots in different energy bands with different peak intensities. 
The peak heights are normalized to be unity. The solid lines represent the average profiles 
superposed for the high intensity shots with the peak count C p > Ct> + 4<r, the crosses for the 
low intensity shots with C p < C^ + Aa, where is the count of steady emission, a=JC p + 
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Figure 3: The superposed shot profiles with the different durations, the solid lines correspond 
to shots with shot durations shorter than 130 ms, and the dash lines correspond to shots with 
durations longer than 130 ms. The peak heights are normalized to be unity. al)-a3) for the 
low to high transition; bl)-b3) for the high state; cl)-c3) for the high to low transition; dl)-d3) 
for the low state. 



is similar too. Therefore, the average shot obtained by superposing all the selected shots in 
our studies can reflect the main characteristics of shots with the different durations and peak 
intensities, especially the rapid variability properties. 

The present method was also applied to the Crab data obtained with PC A/ RXTE. The 
average short time scale ( ~ 10 ms) flare obtained is similar to the one of the two peak average 
pulse profile found by folding with period of 33.4 ms, in duration and rough shape, as shown 
in Figure 4. It can be found some smaller flares which reflect the periodic double-peak pulse 
structure of Crab. It means that the flares selected, using our method, are mostly pulses of 
Crab. But the double-peak structure found by using our method is not as clear as by period 
folding, it is due to most of periodic pulse selected by us has only one peak of its double-peak 
pulse which has significantly positive fluctuation to be selected by our criterion. It shows that 
our method can select the real rapid flares and obtain the rough average shape of flares which 
reflect some properties, such as main flare's rise of the flares in a real light curve obtained by 
PC A/ RXTE in time scale of ~ 10 ms. 



6 



4.6x10 




3.8x1 



1 20 
Time (ms) 



30 




700 



600 EL 



-0.10 



-0.05 0.00 0.05 

Time (s) 



0. 1 



Figure 4: The comparison between the superposed flare profile and the pulse profile in Crab 
data with PCA/ RXTE. a) Pulse profile obtained by period folding, b) Superposed flare profile 
obtained by our method. 
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Figure 5: The superposed shot profiles. The peak heights are normalized to be unity. The 
profiles al-a3 are obtained by superposing 537 shots of the low to high state transition in the 
energy bands of 2-6, 6-13 and 13-60 keV respectively; bl-b3 by 438 shots of the high state in 
the 2-6, 6-13 and 13-60 keV bands respectively; cl-c3 by 268 shots of the high to low state 
transition in the 2-6, 6-13 and 13-60 keV bands respectively; dl-d3 by 564 shots of the low 
state in the 2-5, 5-13 and 13-60 keV bands respectively. 



3 Results 



3.1 Temporal Variability 

The average profiles of X-ray shots have been obtained from PCA / RXTE observations in the 
low-to-high transition, high, high-to-low transition and the low states of Cyg X-l by our im- 
proved superposing procedure described in the last section and shown in Figure 5. In the 
transition states, the front of shot could be simply represented by a single exponential. How- 
ever, in the high and low states, a sum of two exponentials which correspond to the rapid and 
slow variability respectively, is necessary. The formula for fitting is 

— -l>'-C (*<0) (2) 
r+ + B+ e "+ + C + (t > 0) 




8 



where r_ and r + correspond to the rapid process of shot; and r] + correspond to a slow 
process of shot; t is the time relative to the shot peak. The best-fit parameters for the shots, 
within 0.45 s or 450 bins with 1 ms time bin on both sides of shot peak, are summarized in 
Table I. It shows that the rise of the average shot is slower than its decay. And the difference 
between the slopes of rise and decay in the high energy band is more significant than that in 
the low energy band. 

The comparison of shot features in the different states and energy bands are shown in 
Figures 6 and 7. It can be found that the profiles of shots depend upon state and energy. The 
shot duration is the shortest in the high state and the longest in the low state. The shot in 
the high energy band is narrower than that in the lower energy bands. We also found that the 
time constant for the rapid decay in the high energy bands is always about 6 ms in the three 
states. 

The distributions of the time lags between the shots of the high and low energy bands in 
the different states are shown in Figure 8. The time lag is determined by the maximum of 
the cross-correlation function between the shot profiles in the two energy bands. The average 
lags are listed in Table 2, where the errors are estimated by the bootstrap method (Diaconis 
& Efron, 1983). It can be seen from Figure 8 and Table 2 that there exists, on the average, a 
delay of high-energy photons relative to low-energy ones in the X-ray shots of Cyg X-l in each 
state, and in the transition sates the distributions of time lags are most asymmetric and the 
average lags are more larger than that in the high and low states. 

3.2 Spectral Properties 

Figure 9 shows the hardness evolution of average shots in the different states. The zero points 
of time in Fig. 9 correspond to the shot peaks. To study the hardness during a shot, the shot 
peaks in each energy bands have been taken as the bins which are the maximum count bins in 
the summed up shots of every energy bands in 2-60 keV. And the time bin of light curves are 
set to be 10 x 2~ 10 s to suppress statistical fluctuation. From Fig. 9 one can see that in the 
high state the hardness during a shot is higher than that of the steady emission around the 
shot . However, in the low and transition states, the hardness is lower than that of the steady 
emission and the hardness reduces significantly about 100 ms before the shot peak and rises 
rapidly to the maximum near the peak. This result is similar to that obtained by NMK. The 
hardness ratios for the net shots after subtracting the steady emission in the different states 
are shown in Figure 10. It shows that after a drop the hardness rises gradually and reaches to 
a flat top during the shot in the low or transition states. 

We take the region from —0.4 s to —0.2 s and from 0.2 s to 0.4 s from the shot peak as the 
off-shot region and the interval between —0.2 s and 0.2 s around the peak as the shot region. 
The net-shot spectra can be obtained by subtracting the off-shot spectra from the shot spectra. 
The average net-shot spectra in a state can be obtained by merging all the net-shot spectra in 
the state and the average off-shot spectra obtained by merging all the off-shot spectra. Figure 
11 shows the obtained net-shot and off-shot spectra in each state separately. Figure 12 shows 
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Figure 7: Autocorrelations of average shot profiles in the different energy bands. 

the net-shot spectra in the different states for comparison. It can be seen that the spectra of 
net-shot in the low state is the hardest one. 

4 Discussion 

The shot durations found from observation of PCA/RXTE are significantly shorter than that 
obtained with Ginga data by NMK. The longer time constant of shots found from PC A/ RXTE 
observations is similar to the shorter time constant obtained by NMK, while the shorter time 
constant obtained by our studies is shorter than the time resolution of Ginga. The typical 
interval between the two successive shots is about 0.4 — 0.5 s. The longer duration found by 
NMK could be due to the relatively lower time resolution in the Ginga data and the overlapping 
effect of shots in their superposing process. The average shot profile obtained by superposing 
the simulated data shows that such methods tend to smooth the structure of individual shot. 
Although the profiles of shot are not identical as shown in Figure 3, the shots with different 
profiles share some common features, such as the hardness evolution properties and rapid decay 
within ~ 20 ms of shot peak, are reflected by the average shot. 

The shot features vary in different energy bands even in the same state. The duration of 
shot in the high energy band is shorter than that in the low energy bands, as seen in Figure 
6, and it is more significant in the transition and high states than that in the low state. It is 
not consistent with the prediction of the models in which X-ray shots are produced by inverse- 
Compton scattering of flares of soft radiation. Because the higher the energy of photons is, 
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Figure 9: The hardness ratio profiles of average shots in the different states. al-a3 for the 
low-to-high transition; bl-b3 for the high state; cl-c3 for the high-to-low transition; dl-d3 for 
the Low state. 
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Figure 10: The hardness ratio profiles of average net-shots in the different states. The dashed 
lines correspond to the average off-shot hardness. al)-a3) for the low-to-high transition; bl)-b3) 
for the high state; cl)-c3) for the high-to-low transition; dl)-d3) for the low state. 
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Figure 11: The net-shot and off-shot spectra with normalized intensity in the different states. 
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the more the scatterings for the low energy photons are needed (Sunyaev & Titarchuk 1980; 
Hua et al. 1997; Bottcher & Liang 1998). The reirradiation of the high energy photons by 
cool matter on disk could be a explanation of that the duration of shot is larger in the low 
energy bands than that in the high energy bands, but it will cause a lose of the near-perfect 
coherence between the low- and high-energy x-ray photons observed in X-ray emission of Cyg 
X-l (Vaughan & Nowak 1997). The propagation model (Manmoto et al. 1996) assumes that 
an X-ray shot comes from a disturbance added onto the advection-dominated disk at large 
radii, which propagates inward and is reflected at the innermost part of the disk as an outgoing 
shock wave. This model can explain some shot features obtained by NMK. In the disturbance 
propagation model the observed main shot properties, i.e. the shot duration is shorter and 
its peak is later than that in the low energy, could be naturely attributed to that the higher 
temperature the region is, the more inner and smaller it is, if the temperature rises gradually 
from the outer to inner region of the disk. But the propagation model expects the resulting 
light profile to be roughly symmetric with respect to the shot peak, i.e., the rise and decay parts 
should be roughly the same, which is not consistent with the asymmetric structure obtained 
from the observation of RXTE presented in this paper. 

Our results show that some temporal and spectral characteristics of X-ray shots of Cyg X-l 
are remarkably different in the different states. The shot duration is the shortest in the high 
state and the longest in the low state. The evolution of hardness during a shot in the high 
state is significantly different from that in the transition and low states. It could be a new 
criterion to distinguish the high state from the other states of Cyg X-l and interesting to apply 
this criterion to other black hole binaries to judge their states. No proposed shot model has 
predicated such a difference in temporal and spectral properties between the high state and 
other states. The factor of more effective cooling in the disturbance has been used to explain 
the pre-peak drop of hardness during a shot in the transition and low states, but it could not 
explain the hardness during a shot being higher than that of the steady emission in the high 
state. It may imply that the environment of the region where the rapid fluctuation happens is 
changed with the spectral state. The observed hard X-ray emission may mainly be produced 
by Comptonization in a hot advection-dominated accretion flow (ADAF) (Esin et al. 1997) or 
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an accretion disk corona (ADC) (Dove et el. 1997). The joint region between the innermost 
orbit of the disk and the hot corona could be the place where most violent disturbances occur 
and X-ray shots come from. In the low or intermediate state, the hot flow or corona can be seen 
as a sphere surrounding the black hole horizon with an extending accretion disk and the joint 
region is the outer part of the hot corona with a lower temperature of electrons. On the other 
hand, in the high state, the ADAF or ADC is restricted to the corona above the disk (Esin et 
al. 1997) and the joint region is then the inner part of the corona with a higher temperature. 
The geometries of hot corona may explain the difference of hardness evolution in the different 
states. The shot duration may relate to the size of the hot corona. The shot durations obtained 
in the different states implies that the size of corona is smaller in the high state than that in 
the transition states, and it is the largest in the low state, which is also consistent with the 
ADAF or ADC model. 

The time lag of high energy photons relative to low energy ones is a useful probe of the 
dynamics of the hard X-ray production region. Although a time delay of X-ray shots in the 
high energy band is, on the average, observed for each state of Cyg X-l (see Table 2), which is 
consistent with results from the Fourier cross spectral analyses, but the features of the distri- 
butions of time lags are complicated and difficult to be explained with any simple production 
model. A futher sdudy of the time lag distributions is in the process. 

We thank James Lochnner for his help in obtaining and reducing the data, Shuangnan 
Zhang for his valuable comments, and Fanjin Lu and Yong Chen for their worthful discussions. 
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Archive Research Center Online Service, provided by the NASA/Goddard Space Flight Center. 
We also thank the anonymous referee for his very useful comments to improve our study. 
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Table 1. Best fit parameters for the superposed shot in different states 
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Table 2. Average hard X-ray lags of shots in different states. 



State Time Lag (ms) 



Low to High 9.9±0.9 

High 3.9±1.0 

High to Low 9.1±1.5 

Low 1.8±0.7 
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